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We discuss the perturbative running of a s in a model independent way. Our analysis contains 
data on the hadronic r decay and hadronic cross sections in e + e~ annihilation between 5 GeV and 
Mz- We determine the color coefficients and the perturbative /3-function of the strong interaction. 
The results are in agreement with QCD and rule out the QCD+light gluino scenarios on the 70.8 
- 93.0% CL. We combine our method with that of the multi-jet analysis at LEP. The combined 
result rules out light gluinos at least on the 99.76% CL, provided that nonperturbative effects are 
not large. 

Asymptotic freedom is one of the most interesting predictions of QCD. In order to study the running of a s , one 
can collect its values at different scales (e.g. r decay, DIS, T decay, e + e _ ) and compare them |IJ]. Another way is 
to analyse a single experiment (e.g. pp, ep or r decay 0]). Nice agreement has been found between experiments and 
QCD. Some differences between the results obtained from low-energy and high-energy experiments led to a number 
of speculations; e.g. that the apparent slower running of a s could be due to additional light fermions of the theory. 

Supersymmetric phenomenology deals normally with sparticles of masses O(l00GeV). The only exception is the 
light gluino || with mass < 1.5 GeV (window I) §,§ and 3 - 5 GeV (window III) [§. We will discuss window I 
and window III separately. Since such a light particle influences the running of a s the comparison of high and low 
energy experiments could open or close the light gluino window [Q— ||. A consistent analysis must contain the virtual 
gluino effects not only in the running (technically in the /3-function), but in all loop diagrams, extracting a s from the 
experimental results, too ||. 

One can even determine the /3-function of the strong interaction as done in ref . P] . In order to extract the three-loop 
coefficients one needed large values of a s , that is the small energies of r decay. The result is in good agreement with 
QCD. 

It has been suggested Jl0| to look for gluinos in four-jet events. The LEP collaborations determined the color 
coefficients (Ca/Cf and Tf/Cf) in multi-jet events. The results (c.f. |j]]) are in good agreement with QCD. An 
extension of QCD including light gluinos would result in a somewhat different Tf/Cf- The four-jet cross-section is 
known only at tree-level. Without the knowledge of the loop corrections this method alone can hardly give a clear 
answer to the question of light gluinos. 

In this letter we determine the color coefficients and the perturbative /3-function of the strong interaction from 
experimental results. We consider gauge theories with fermions in the fundamental representation and their extended 
versions with light gluinos. We perform an analysis on the three-loop level in the MS scheme. The experimental 
inputs are R T and Rh for energies between 5 GeV and Mz- We combine our method with that of the multi-jet analysis 
at LEP. 

One of the clearest ways to study the running of the coupling, thus the /3-function, of the theory would be to study 
one given experimental quantity at different energies. We chose Rh's at energies larger than 5 GeV, which are very clear 
quantites both experimentally and theoretically, due to their minimal nonperturbative corrections. Moreover, these 
are the only quantities known to O(a^) relative to the leading order, thus only a minimal scale ambiguity is present. 
In order to study the running of a s one needs precise Rh measurements for a large energy region. Unfortunately, Rh 
measurements are limited by their small statistics for energy scales below Mz- For this reason we have also considered 
another, strongly related (still not very low scale) quantity R T in our analysis. Since we are aware of the theoretical 
criticism on a s determinations based on R T we include the effects of all known uncertainties in case of QCD 0- 
p3[ . The nonperturbative estimates to R T in the light gluino case might be unreliable. Therefore, we do not include 
R T in our analysis for window I gluinos. (Window III gluinos are much too heavy to contribute to R T .) 

Let us suppose that the strong interaction is described by a gauge theory based on a simple Lie-group, which fixes 
the color coefficients of the theory (Cf,Ca,Tf)- The three basic processes in the theory are: gluon bremsstrahlung 
from a quark, splitting of a gluon into two gluons and splitting of a gluon into two quarks. To lowest order their 
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amplitudes are proportional to a universal coupling and to Cp, Ca and Tp, respectively. E.g. for an SU(N) gauge 
group one has Cf — (N 2 — l)/(2iV), Ca = N and Tp = 1/2. The /3-function of the strong interaction, hadronic 
cross-sections and widths are calculated in terms of these coefficients and the active number of fermions. We will look 
for a set of parameters which describes the experimental data most accurately and give the corresponding confidence 
level (CL) regions. The outcome -that is the best fit- does not necessarily predict a meaningful theory. Nevertheless, 
it could tell the difference between QCD and QCD+light gluino scenario. 

The /3-function: da(n) / d\og fi — —(3oa 2 — /3ia 3 — /^a 4 is known upto three-loop order. We factorize out the 
inconvenient 7r-s by using a = Cpa 8 /(2Tt). Introducing x — Ca/Cf and y = Tf/C'f gives /3q = llx/3 — 4(n/y + 
xrig/2)/3, where n/ is the number of active flavours and rig is the number of active gluinos. The expressions for /3i 
and are more complicated [|lj,^5|. In the running of a s we follow ref. for threshold-effects. We solve the above 
renormalization group equations exactly. 

The other important quantity in the analysis -the hadronic cross section in e + e~ annihilation via a virtual photon- 
is known on the three-loop level, too. 

R 1 = a(e + e~ — > 7 — > hadr) / a {e + e~ — > 7 — » /z/x) (1) 
= 3E/<#(1 + Kia + K 2 a 2 + K 3 a? + ...), 

where K\ — 3/2. K%, K3 is more complicated and also known for arbitrary Cf, Ca and Tp Jl4],[l^]. The hadronic 
decay ratio of the r is defined: R T = r(r~ — > v T + hadr.) /T(t~ — > v T e~ D e ). Its perturbative value is strongly related 
to R 1 . For a recent analysis and a review on R T see eg. and pl| . The hadronic decay width of the Z boson 
Rz = T(Z — » hadr.) /T(Z — > is again a function of R 1 and known at the three-loop level. At the energy scales 

of the analysis important electroweak and mixed corrections appear. We include these corrections. 

The experimental values for R T are given in jl9| with an average of R r — 3.616 ± .02. Averaging for the four LEP 
experiments and three leptons p0] gives Rz — 20.778 ± .029. The references for hadronic cross sections at energies 
below Mz are taken from We have collected all the existing published data, and some unpublished results, too. 
Some of them were binned by the experimental groups. The total number of data points included in our analysis is 
182. 

Combining results from different experiments is a delicate question, with major problems, i) The use of the 
radiative corrections is not unique, reflecting the state of art at the time of publication, ii) The results do depend on 
the mass of the Z boson and top-quark -which was assumed to be approximately 20 GeV at the early eighties, iii) 
The measurement of the total cross section was not performed in the full phase space, and Monte-Carlo acceptance 
calculations based on different assumptions were done. We corrected for the first two problems, correction for the 
third one is practically impossible. For further details see [p2|- |23| ]. 

There are different sources of uncertainties in the determination of a s . We treat them in a unified manner. We 
add the systematic errors linearly. The total systematic error estimates and the statistical errors are combined 
quadratically. The overall normalization errors within one experiment are correlated. We include these correlations. 
We minimize x 2 = A T V~ 1 A, where A is an n- vector of the residuals of Ri — Rfn for the n individual results and V 
is an n x n error matrix. In V, the diagonal elements Vu are the squares of the total errors for the ith measurement 
and the off diagonal elements Vij correspond to the correlations between the ith and jth measurements. For different 
points from the same experiment Vij is given by the product of the normalization errors p^] . The separation of 
the systematics into point-to-point and overall normalization error is ambiguous. Our results are rather insensitive 
to a moderate change of this splitting in a given experiment. For some experiments the above separation was not 
even explicitly given. We checked that our results are stable against even large variations of the splitting in these 
cases. However, the unreasonable extreme case of totally uncorrelated fitting would give rather different results with 
much more predictive power. We have assumed that the results of different experiments are uncorrelated. In order 
to perform a further consistency check we included an additional hypothetical overall error of 1% in the correlation 
matrix. The change of the result turned out to be negligible. Since the experimental groups usually give their results 
with binning in fixed energy intervals we checked that further binning had practically no influence on our results. 
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FIG. 1. The 68.3% CL regions for x = Ca/Cf and y = Tf/Cf- The QCD value is represented by a square. 

Due to the relatively poor statistics for the e + e~ experiments below Mz the theoretical uncertainties are dominant 
only for R T and non- negligible fori?^. Estimating the nonperturbative and perturbative errors in T-decay we use 
the results of [0,01. For a fixed order QCD calculation we assume that the error is equal to the last computed 
term. Including all the errors we get R T = 3.616 ± .143, which corresponds to a s (M T ) = .335 ± .053 and Rz = 
20.778 ± .0387, which corresponds to a s (Mz) — .123 ± .006. We assume that currently incalculable and/or model 
dependent nonperturbative corrections are negligible or correctly estimated. If the nonperturbative corrections turn 
out to be larger than our estimates, the confidence levels for light gluino exclusion to be determined below should be 
lowered. Note that both error estimates are certainly very conservative. 

Having included all the errors one can determine the best fit values and the CL regions for x = Ca/Cf and 
y = Tf/Cf with or without light gluinos. 

Fig. 1 contains the 68.3% CL regions for theories without gluinos and with window III gluinos. As it can be seen 
the experimental results are in agreement with QCD; however, they rule out the QCD+light gluino scenario on the 
93 (90.7)% CL's for Mg = 3 (5) GeV, respectively. The variation of the fixed x 2 boundaries is quite small for the 
above mass region. For binning level 50 (d.o.f. 48) the \ 2 values are 37.37 without and 37.94 with gluinos (Mg = 3 
GeV) for the best fits. As mentioned before, for window I gluinos we do not include R T in our study. In this case we 
have a 70.8% exclusion only. 
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FIG. 2. The 68.3% and 95.4% CL regions for (3 and fa. 

Note the important difference between the presentations of our result and that of the 4-jet analysis for the case 
with light gluinos B. In our case the best fits and the CL contours for the x, y variables are determined separately 
for the theories with and without light gluinos. These results are compared for both theories with the QCD values: 
x = 2.25, y — .375. Due to the simplicity of the tree-level treatment of the 4-jet analysis in ref. M, it was possible 
to produce the same best fits and contours for both theories and compare this unique result with x, y above in the 
theory without light gluinos, and with x and with an effective y e ff = .6 for the gluino extended theory. The 4-jet 
results can easily be expressed in terms of our variables. 

Two remarks are in order: i) the 4-jet analysis is based on tree-level calculations, whereas our method contains 
corrections up to three-loops. (In jlj also the one-loop analysis to 2- and 3-jets is done; however, the errors on x, y 
are huge.) ii) our method contains the running of a s , too, thus it gives a simultaneous check of the /3-function. 

Assuming [ pij that the underlying gauge group is SU(3), we may consider the number of gluinos as a free parameter. 
The best fits give n g = —.638 ± 1.17 for window I and n g = .0078 ± .52 for window III. The CL's are given by Bayes' 
theorem allowing only non-negative integer n g : 72.5% for window I and 87.7% for window III. (For modeling data 
with bounded physical region: n g > 0, and for application of Bayes' theorem see e.g. Sec. 28 of |pl| .) Finally, one 
can also fix the number of gluinos and determine a s (Mz)- For rig = (i.e . QCD) we get .124 ± .004, for n g = 1 we 
get .129 ± .006 (window I) and .132 ± .006 (window III). As expected [£lj the QCD value is slightly larger than the 
world average. 
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FIG. 3. The 68.3%, 95.4% and 99.7% CL regions of the combined analysis for the theory with window III light gluino 
(Mg — 3 GeV) and the 68.3% CL region for the theory without gluino. 

Fig. 2 shows the /3-function coefficients with 68.3% and 95.4% CL regions for models without gluinos at Mz- These 
curves represent merely a transformation of the curves in the {x, y) plane, they result in similar conclusions for QCD. 
As it can be seen the value of (3o is quite well constrained; however, not even the sign of (3i could be predicted. This 
result reflects the fact that in the energy region studied the running of the coupling is determined almost exclusively 
by the one-loop term, proportional to [3q. An interesting possibility is to fix the cross sections, /3i, /?2 and the flavour 
dependence of /3q to their perturbative QCD values and extract /3o from the fits. The result is /?o = 5.70 ± .64 to be 
compared with the theoretical value of 5.75. Following (|] we fixed /?o and /3i to their theoretical (QCD) value and 
performed a fit with 02 as a free parameter. The allowed region for @2 contains the theoretical prediction, but the 
error is an order of magnitude larger than the theoretical value. 

Comparing our Fig. 1 and Fig. 10 of § -taking proper care of using the same variables- one observes that the 
overlap of the 68.3% CL regions for the two analyses is quite small. This opens the possibility for a much stronger 
restriction on the theories without and with gluino. We have parametrized the % 2 of the multi-jet analysis as a 
function of x and y, and performed a new fit adding it to our y 2 function. In the actual analysis we have included 
the new results in (2(|. Since the 4-jet analysis is based on a tree- level calculation we have included the unknown 
higher order QCD and mass effects. We assumed much larger uncalculated higher order corrections than used in the 
experimental papers. Due to them we enlarged the axes of the error ellipse by 12% of the theoretical x and y values 
(relative correction of 0(a s ).) For the 4-jet analyses we have used the mass effects calculated by pq| . The results of 
the combined analysis for QCD and window III gluinos are presented in Fig. 3. Our result is consistent with QCD + 
no gluino scenario, while the theory with window III light gluino is excluded on the 99.99 (99.89)% CL for Mg = 3 (5) 



GeV. The window I gluino is excluded on the 99.97% CL. 

Again one can fix the underlying gauge group to SU(3) and consider the number of gluinos as a parameter to be 
fitted. We get n~ g = -.35 ± .33 for window I and n g = -.156 ± .27 (n§ = -.197 ± .32) for window III M s = 3 (5) GeV 
gluinos. Using Bayes' theorem assuming that n g is non-negative integer number, the exclusion CL's are 99.96% for 
window I and 99.96 (99.76)% for window III gluinos. 

Leaving out R T from the window III combined analyses, the CL values change only negligibly. 

In summary, we have presented a method to analyse the running of a s from M T to Mz- We have determined the 
group coefficients of the theory, which are in agreement with QCD. The /3-function coefficients are also given with 
a similar result. Furthermore, our results rule out the QCD+light gluino scenario on the 70.8% and 90.7-93.0% CL 
for window I and III, respectively. Combining the analysis of the present letter and the method using multi-jet data 
(cf. ref. [|l 26 1, however treating the theoretical errors more conservatively than the experimental groups) we have 
got an order of magnitude smaller CL regions for the allowed group coefficients, which supports QCD within la and 
excludes the QCD+light gluino scenario with a 99.97% for window I and at least 99.89% CL for window III gluinos. 
Performing a one parameter fit to determine the number of gluinos, we get n g = —.35 ± .33 and n§ = —.197 ± .32 for 
window I and III, respectively. In terms of CL's this corresponds to 99.96% and 99.76%. Needless to say, our results 
are valid assuming the light gluino extension of the standard model. We have not investigated whether additional 
new physics could change our exlusion CL limits. 

The details of the above analysis, the convergence features and the statistical analysis will be presented in a 
forthcoming publication |27f] . 
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